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SIMULATION  OF  DISLOCATION  AND  TRANSFORMATION  PLASTICITY  IN  SHAPE 

MEMORY  ALLOYS 

Summary  of  results 


There  has  been  developed  a  model  of  deformation  of  shape  memory  alloys  based  on 
the  methodology  of  the  background  microstructural  approach  proposed  by  Prof. 
V.A.Likhachev  and  named  the  structure-analytical  theory  of  strength.  Two  scales  of  plastic 
deformation  are  modeled.  The  first  is  the  micro-plastic  deformation  caused  by  the  interphase 
stresses  arising  at  the  accommodation  of  martensite;  the  second  is  the  active  plastic 
deformation,  produced  by  external  mechanical  loading.  Constitutive  equations  for  alloys  of 
T1N1-  and  FeMnSi-types  have  been  formulated  and  computer  programs  composed.  The  model 
can  describe  such  phenomena  as  the  incompleteness  of  the  strain  recovery  at  thermocycles, 
plasric  deformation  of  austenite,  repeated  two-way  shape  memory  effect,  produced  both  by 
micro-plastic  and  macro-plastic  deformation.  Calculated  dependencies  are  in  a  good 
qualitative  agreement  with  the  available  experimental  data. 
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Introduction. 

At  present  the  problem  of  modelling  of  the  mechanical  behaviour  of  shape  memory 
a  oys  (SMAs)  becomes  more  and  more  significant  because  many  new  applications  in  various 
fieWs  of  technology  and  medicine  [1,3]  are  being  developed.  A  computer  simulation  helps  to 
check  the  workability  of  a  device,  optimize  its  parameters.  In  1980s  -  1990s  there  were 
^l0Plmany  theoreticaI  approaches  which  can  describe  the  mechanical  behaviour  of 
SMAs.  They  may  be  divided  into  two  groups.  The  first  consists  of  phenomenological 
macroscopic  theories  [4-8]  which  establish  a  relation  between  stress,  strain,  internal 


parameters  and  their  rates  mainly  on  the  grounds  of  observations  of  their  mechanical 
behaviour.  The  advantages  of  such  theories  are  their  compactness,  a  comparatively  easy 
procedure  of  determining  of  the  material  constants,  relatively  small  amounts  of  computation. 
The  shortcoming  is  their  inter-  or  extra-  polation  character,  less  predictive  ability.  The  other 
group  is  formed  by  microstructural  theories  [9-12].  They  are  based  on  crystallographic  and 
physical  data  about  the  structure  of  the  material  and  mechanisms  of  deformation  and  take  into 
account  geometrical  and  physical  peculiarities  of  deformation  mechanisms.  This  allows  to 
develop  rheological  models  invariant  of  the  external  loading,  to  make  predictions  about  the 
behaviour  of  a  particular  alloy  under  the  action  of  a  complex  thermomechanical  load  most 
accurately  and  physically  grounded.  At  the  same  time  in  such  theories  it  is  more  difficult  to 
calibrate  material  constants,  big  computation  resources  are  required  and  this  impedes  their 
use  for  .solving  boundary-value  problems.  Still  the  necessity  to  describe  different  regimes  of 
an  SMA  application  in  the  frames  of  one  general  physically  justified  approach  makes  the  use 
of  micro-structural  theories  very  attractive.  In  the  present  project  a  background 
microstructural  approach  developed  by  V.A.Likhachev  and  named  "structure-analytical 
theory  of  strength  is  used  as  the  guideline  and  an  attempt  is  made  of  its  expansion  to 
loadings  causing  irreversible  plastic  deformations. 

Basic  constitutive  equations  of  structure-analytical  theory  of  strength 

The  structure-analytical  theory  of  strength  [9,10]  assumes  that  the  main  characteristic  of  a 
crystal,  determining  many  of  its  properties  and  particularly  deformation  mechanisms,  is  the 
crystallographic  orientation  to  which  the  orientation  of  slip  systems  is  related.  A  polyciystal 
consists  of  grains  differing  by  such  orientations,  therefore  to  obtain  macro  deformation  it  is 
reasonable  to  use  orientation  averaging  of  deformations  of  grains.  This  averaging  may  be 
supplemented  with  the  statistical  one  on  such  parameters  as  the  width  of  thermomechanical 
hysteresis.  The  averaging  procedure  is  organized  in  a  way  to  reflect  the  grain  structure  of  the 
material  [13].  In  each  of  1,2 grains  of  the  representative  volume  of  a 
polycrystalline  body  the  deformation  related  to  the  crystallographic  basis  of  that  grain  is 
characterized  by  a  matrix  fir( a)  .  The  macroscopic  phase  strain  is  obtained  by 
reducing  matrices  fir(a>)  to  the  laboratory  basis  and  averaging  on  all  grains : 

=  in<o)RaeGr{co)R-',  (1) 

where  Ra  is  the  rotation  matrix  transforming  the  laboratory  basis  into  the  basis  of  grain  co 

ancj  f(°)  1S  volume  fraction  of  this  grain.  The  volume  fraction  <jM  of  martensite  in  a 
grain  is  calculated  by  the  same  averaging  procedure: 

0M  =  Z  f(o))0Gr(co ),  (2) 

On  the  total  in  each  grain  there  can  exist  1,. N  variants  of  martensite.  A  variant  is 
characterised  by  quantity  ®n,  such  that  (1/N)®n  is  its  volume  fraction  in  the  grain,  and 
lattice  deformation  Dn  by  which  it  was  obtained  from  the  parent  phase.  If  the  matrix  D  for 
one  of  the  variants  is  known,  the  corresponding  matrices  Dn  for  other  variants  can  be  found 
by  rotating  D  by  operators  Pn  belonging  to  the  cube  symmetry  group: 

At  =  PnDPn~l,  n  =  1,...,  N. 


(3) 


In  each  grain  the  average  phase  strain  accumulated  due  to  lattice  deformation  in  the 
course  of  the  transformation  related  to  the  crystallographic  basis  and  the  volume  fraction  of 
martensite  are  equal  to 


1  N 

*Gr=772X’ 


N 


n-  1 


(4) 

(5) 


For  simplicity  ®Pr,  ®n,  eph,  etc.  here  and  further  are  written  instead  of  <x£r(co),  &n(co), 
^(co),  etc. 

To  formulate  the  law  of  variation  of  &n  we  consider  the  thermodynamic  equilibrium  of  the 
two-phase  austenite-martensite  system.  We  introduce  the  thermodynamic  force  Fn{  which  is 
the  difference  of  specific  Gibbs  energies  of  the  n-th  variant  of  martensite  and  austenite.  An 
approximate  expression  for  Fnl  is 


Fnf  =  (Tn*~  T0)iq/r0) ,  with  Tn*  =  T+  (T/q^vDn,  (6) 

where  Tn  is  the  effective  temperature,  T0  is  the  temperature  of  the  equilibrium  of  the  two 
phases  when  there  is  no  stress,  qo  is  the  latent  heat  (qo< 0),  t  is  the  stress  and  double  dot 
sign  denotes  the  double  contraction  of  tensors.  It  is  generally  assumed  that  a  moving  interface 
overcomes  a  resistance,  due  to  the  energy  of  martensite  crystal  nucleation  and  other  barriers. 
Here  we  call  it  friction  force  F  fr  and  assume  that  it  has  a  constant  absolute  value  and 
opposes  the  movement  of  a  phase  boundary.  Besides  friction,  there  appears  a  counter-force 
Fn  due  to  the  interaction  between  martensite  crystals  and  austenite  matrix.  We  postulate  as 
the  first  approximation  that  this  force  can  be  calculated  as: 


Fn~Fm^n  (7) 

where  /um  is  a  constant. 

At  the  thermodynamic  equilibrium  there  is  the  balance  of  forces: 

Fnf  -  Fn  ±  ,  (8) 

where  plus  sign  must  be  taken  for  the  direct  and  minus  —  for  the  reverse  transformation. 
Constants  Ffr  and  /um  can  be  determined  from  the  kinetics  of  the  martensitic 
transformation  when  no  stress  is  applied  and  they  are:  Ffr  =  (~qo)'(Af  -  Ms)/(Af  +  Ms), 

Fm  ~  (~qO>2(Ms  ~MJ)/(Af+  Ms)  ,  Mf ,  Ms ,  As ,  Af  being  as  usual  the  characteristic 
temperatures  of  the  transformation.  From  (7)  and  (8)  differential  d®n  can  easily  be  found: 

d®n=  dFn1/ Fm  •  (9) 

Initial  conditions  are  based  on  the  fact  that  at  temperatures  above  yfyand  in  stress-free 
state  no  micro  deformations  have  begun  and  for  all  martensite  variants  we  have  ®n~  0 


Theoretical  studies  of  different  variants  of  the  structure-analytical  model  [11,13]  has 
shown  the  capability  of  this  approach  to  correctly  describe  the  basic  features  of  SMA 
behaviour  at  different  regimes  of  thermal  and  mechanical  loading. 

Modelling  of  martensite  accomodation  effect  on  mechanical  behavior  of  TiNi  type  shape 

memory  alloys  [3*] 


The  first  step  to  expand  the  field  of  applicability  of  the  model  was  to  take  into  account 
micro  plastic  deformation  and  stress  produced  by  the  accommodation  of  martensite  [14] 
which  plays  a  big  role  in  the  pecularities  of  SMA  deformation.  In  particular,  it  has  been 
shown  that  accomodation  processes  are  connected  with  such  phenomena  as  the  decrease  of 
the  phase  yield  limit  at  cyclic  loading  [15].  They  are  also  responsible  for  such  phenomena  as 
incomplete  recovery  of  strain,  strain  accumulation  at  thermocycling  and  repeated  two-way 
shape  memory  effect. 

To  describe  the  deformation  produced  by  the  accommodation  of  martensite  consider 
equation  (4)  as  a  decomposition  of  the  strain  £&  relative  to  some  generalized  basis  {Dn} 
with  components  0n  and  which  we  therefore  shall  interprete  as  the  measures  of  strain.  As 
the  micro  plastic  strain  e  mP  caused  by  the  accommodation  of  martensite  is  connected  with 

the  same  martensite  crystals  as  the  phase  strain,  we  postulate  that  s  ™P  can  be  expressed  in 
the  form  similar  to  (4). 


i  N 

=  —  y  0"d 
n  h  "  n 


(10) 


where  0yP  is  the  measure  of  the  micro  plastic  strain  produced  at  the  accommodation  of 
crystals,  belonging  to  the  jj-th  variant  of  martensite  in  a  given  grain. 

As  the  accommodation  leads  to  a  partial  relaxation  of  interphase  stresses,  we  assume  that 
the  resistance  force  Fn  decreases  by  the  quantity  am  0J>  (in  [14]  it  was  taken  pm=am  ). 

So  we  must  change  equation  (7)  for  resistance  force  to 


F n  Pm<&n~clmFrP  (11) 

We  note  that  after  (11)  has  been  accepted  instead  of  (7),  formula  /um  given  above  is  no 

more  valid  as  the  kinetics  of  the  transformation  is  different  when  martensite  accomodation 
takes  place. 

The  evolution  equation  for  0yP  is  obtained  from  a  hypothesis  that  differentials  d0yP 
and  dFn  are  related  in  the  same  way  as  the  differentials  of  strain  and  stress  in  the  classical 
flow  theory  with  linear  hardening.  Then  one  can  easily  show  that  when  the  force  Fn  exceeds 
some  yield  limit  Frf  ( 


d0yP  -  k  d0n  H( |  Fn |  -FfjF)  H(dFns\gnFn),  (12) 

where  H  is  the  function  of  Heaviside,  k  is  the  coefficient  characterising  hardening.  If  k  =  0 
there  is  no  micro  plasticity,  if  k  -»  1  there  is  ideal  micro  plasticity  without  hardening.  As  the 
counter-force  Fn  grows  (decreases)  only  at  the  growth  (decrease)  of  martensite,  which  in  its 


turn  can  occur  only  at  the  growth  (decay)  of  the  thermodynamic  force  Fnl,  we  have 
sign (dFn)  =  sign (dFyf)  and  H(dFnsignFn)=  H(dF„t  signFw).Then 

d0n=  dFyf/  {Mm  [1‘  I  Fn  |  -F^)H{dFy{  signF„)]}  (13) 

When  micro  plastic  deformation  is  in  progress,  the  yield  limit  grows  at  the  same  rate  as  the 
force  Fn{  whose  increment  is  equal  to  that  of  the  thermodynamic  force. 

Initial  conditions  are  based  on  the  fact  that  at  temperatures  above  ^yand  in  stress-free 
state  no  micro  deformations  have  begun  and  for  all  martensite  variants  we  have  F^  =  F(jy  = 
const ,  =  0n=  0 . 

When  micro  plasticity  is  taken  into  account,  the  strain  in  a  grain  is  calculated  as 

-L£  D,(<P. +  <!>;),  (12) 

-/V  /i=i 

and  the  total  macro  strain  is  delivered  by  averaging  (1).  For  numerical  calculations  the 
relations  given  above  can  be  reduced  to  a  system  of  QN  ordinary  differential  equations. 

Examples  of  modelling 

Deformation  simulation  was  carried  out  for  an  object  consisting  of  120  grains  of  different 
pseudo-random  orientations.  The  following  values  of  material  constants  have  been  used:  Mf 
=  280  K,  Ms  =  300  K,  As  =  340  K,  Af=  360  K,  T0  =  330  K,  q0  =  -150  MJ/m3.  The  value 
F0y  ~  13  MPa  has  been  chosen  in  such  a  way  that  microplasticity  could  take  place  at  cooling 
of  specimen  stressed  by  30  MPa.  Parameter  k  has  been  varied.  Calculations  have  been 
performed  for  cases  of  pure  shear  and  uniaxial  tension.  Heat  expansion  has  not  been  taken 
into  account.  Microplastic  deformation  at  accommodation  of  martensite  strongly  influences 
mechanical  behavior  of  shape  memory  alloys  and  especially  at  cyclic  variations  of 
temperature  or  stress.  Because  of  this  we  have  chosen  to  model  strain  accumulation  and 
recovery  by  a  specimen  loaded  with  a  constant  stress  of  100  MPa  and  subjected  to 
temperature  cycles  through  the  full  interval  of  the  phase  transformation.  Fig.  1  presents  the 
kinetics  of  this  process  for  k  =  0.4.  Maximum  strain  saturates,  strain  unrecovered  in  a  cycle 
decreases  and  the  hysteresis  loop  stabilizes  in  approximately  20  cycles  (Fig.2). 

■50 


Fig.l  Strain  vs.  temperature  dependence 
at  thermocycling  in  the  full  interval  of 
the  martensitic  transformation  under  the 
stress  of  100  MPa.  Numerals  on  the  plot 
indicate  numbers  of  cycles. 
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Fig.2.  Calculated  strain  at  a  temperature  below  Mf  (a)  and  the  unrecovery  of  the  strain  (b)  in 
an  N-th  cycle.  Thermal  cycling  was  carried  out  through  the  full  interval  of  the  martensitic 
transformation  under  the  stress  of  100  MPa  symbols  in  fig.2b  correspond  to  the  values  of  k, 
given  in  fig.2a. 

An  important  place  in  the  studies  of  the  work  ability  of  SMAs  is  occupied  by  the  following 
thermomechanical  cycle:  cooling  of  a  sample  from  austenite  to  martensite  under  some  load  P, 
application  of  an  additional  load  Q,  heating  to  a  temperature  higher  than  A/ and  removing  the 
additional  load  Q  to  have  the  initial  load  P.  Figures  3,  4  show  the  result  of  such  numerical 
experiment  when  P  =  100  MPa,  Q  =  200  MPa.  Stabilization  of  the  loop  in  this  case  is  also 
achieved  in  about  20  cycles.  Strain  unrecovered  in  a  cycle  becomes  small  and  almost 
constant.  The  maximum  strain  accumulated  up  to  the  A-th  cycle  slowly  grows,  these  results 
are  in  a  good  agreement  with  the  data  of  the  direct  experiment  [21].  An  application  of  an 
additional  load  Q  at  a  low  temperature  causes  no  unelastic  strain  as  the  stress  produced  by 
P+Q  ^  insufficient  for  the  reorientation  of  martensite.  As  the  temperature  rises,  the 
reorientation  limit  stress  decreases.  At  some  temperature  lower  that  Afihe,  reorientation  of 
martensite  starts.  Then  a  typical  “hump”  appears  on  the  strain  recovery  curves  in  fig.3. 


Fig.  3.  Calculated  strain  versus 
temperature  diagram  at  thermocycles 
with  cooling  under  the  stress  100 
MPa  and  heating  under  300  MPa 
(torsion).  Numerals  indicate  numbers 
of  cycles. 


thermocycles  depicted  in  fig.  3. 

One  of  the  properties  of  SMAs  which  until  now  could  not  be  satisfactorily  described  in  the 
scheme  of  the  structure-analytical  model  is  the  repeated  two-way  shape  memory  formed  by 
cooling  under  a  stress  or  by  active  straining  in  martensite.. 


Fig.  5.  Calculated  strain  versus  temperature  Fig.  6.  Strain  versus  temperature  diagram  of 
diagram  at  cooling  and  heating  under  100  MPa  active  deformation  of  a  sample  in  martensite 
(curve  1),  unloading  and  a  stress-free  state  and  unloading  (curve  1),  heating  (curve 

thermocycle  (curve  2).  2)  and  subsequent  stress-free  thermocycle 

(curve3). 

The  mechanism  of  this  phenomena  may  lie  in  the  appearance  of  micro  stress  fields  as  the 
result  of  the  incompatible  plastic  strain  accompanying  the  growth  of  martensite  crystals 
Figures  5  and  6  present  the  results  of  modelling  of  this  effect.  The  model  demonstrates  a 
reversible  strain  variation  which  corresponds  to  experimental  data  [16].  According  to  existing 
notions  [17]  plastic  deformation  at  martensite  accommodation  plays  an  important  role  in 
training  of  the  material  by  thermocycling  under  a  load,  fatigue  damages,  the  formation  of 
eigenstresses  responsible  for  the  adaptation  of  the  material  and  formation  of  the  repeated  two- 
way  shape  memory.  The  last  thesis  is  confirmed  by  the  fact  that  the  repeated  shape  memory 


appears  only  after  a  strong  enough  straining  of  the  object  by  a  load  exceeding  a  critical  level 
when  there  is  no  complete  strain  recovery.  At  the  same  time  all  these  considerations  were 
mostly  pure  qualitative  and  did  not  allow  a  passage  to  engineering  calculations.  The  proposed 
model  sets  a  formal  ground  for  such  notions  and  the  agreement  of  calculated  results  with 
experimental  ones  confirms  the  validity  of  this  approach.  As  the  result  it  becomes  possible  to 
model  thermocyclic  training  and  repeated  two-way  shape  memory.  Calibration  of  the  material 
constants  allows  to  obtain  valid  results  for  specific  alloys. 

Modelling  of  shape  memoty  alloys  with  fcc-hcp  transformation  (of  FeMnSi  type)  [1*] 

One  of  the  specific  features  of  fcc^hcp  transformations  is  the  multi-variance  of  both 
direct  and  reverse  transformations  [16,18].  The  result  of  it  is  that  specimens  loaded  by  a 
constant  stress  experience  deformation  when  they  undergo  either  of  the  transformations. 
Thus,  considerable  strain  can  be  accumulated  at  thermal  cycles.  Large  part  of  this  strain  is 
reversible  and  can  be  recovered  at  further  cycles  after  the  stress  is  removed  [19].  This  fact 
suggests  that  the  strain  is  due  to  transformations  rather  than  to  dislocation  movement.  Strain 
accumulation  at  thermocycles  can  be  described  by  the  structure-analytical  theory  of  strength 
which  is  distinguished  by  two  ideas:  (1)  the  existence  of  limits  of  the  size  of  growing 
martensite  crystals  and  (2)  the  possibility  that  the  deformation  of  a  martensite  crystal  at  the 
reverse  transformation  is  not,  in  general,  opposite  to  the  deformation  which  this  crystal  has 
undergone  at  the  direct  transformation.  The  model  based  on  only  these  two  features  predicts 
absolutely  symmetrical  behaviour  of  a  loaded  sample  in  heating  and  in  cooling.  An  attempt  to 
account  for  an  asymmetry  has  been  made  in  [20],  However,  approaches  [10,20]  did  not  take 
into  account  microplastic  strain  arising  at  the  growth  of  martensite  crystals  and  leading  to  a 
partial  accommodation  of  martensite  and  relaxation  of  interphase  stresses.  Because  of  this 
drawback  the  stabilization  of  the  deformation  of  the  model  material  occurs  in  3  -  5  cycles 
rather  than  in  several  dozens  as  in  an  actual  experiment.  Strain  recovery  at  thermal  cycles 
without  stress  can  not  be  described  either.  The  present  work  accounts  both  for  symmetrical 
features  of  fee  <->  hep  transformations  and  microplastic  deformation. 

At  fcc^hcP-transformation  one  of  { 1 1 1 }  fcc  planes  becomes  (0001 )hcp  by  one  of  three 

shears  1/6  <  1 12  >  fcc  on  each  second  plane  making  up  N=  12  variants  of  transformation 

strain  tensor  Dn  and  martensite  quantities  <E>n  .  At  the  reverse  transformation  each  of  three 
shears 

1/3  <  1 120  >hcp  restores  the  initial  orientation  of  austenite.  Thus  all  variants  can  be  divided 
into  four  triplets  (zones)  with  parameters 

^z-r  >  z=\,2,3,4,  (3) 

Jn=3z-2 

characterising  the  amount  of  martensite  belonging  to  a  zone,  (1/4)0Z  being  the  volume 
fraction  of  martensite  of  this  zone.  Parameters  0n  themselves  do  not  have  physical  meaning 
of  martensite  amounts  as  deformations  Dn  belonging  to  one  zone  produce  martensite  of  the 
same  orientation.  Still  &n  play  the  role  of  measures  of  the  phase  deformation.  Total 
deformation  of  a  grain  caused  by  transformation  and  accommodation  of  martensite  is 
presented  by  the  same  formula  (12).  Evolution  equations  for  &n  is  obtained  by  the  same 


considerations  as  for  TiNi  type  alloys.  The  difference  is  that  two  extra  conditions  of 
transformation  finish  exist:  direct  transformation  can  not  occur  to  make  the  volume  fraction 
of  martensite  in  a  grain  OS7"  >1;  martensite  fraction  in  a  zone  (1/4)  ®z  can  not  become  less 
than  zero.  Summarizing  all  this  we  arrive  at  a  final  form  of  the  evolution  equation: 

d0 .  -  {%iHh(|F.  I  -  F„y  )H(-dFn*  signFa )  +  h(f.’  -  |F.  fe  + 

+  ,Gfl  th"  H(lFn  I  -  Fn  )H(rfF>/(gwFn)dFD  }h(i  -  QP)n(0t) ,  (11) 

For  numerical  simulations  the  initial  value  of  the  yield  limit  (F/=10  MPa)  was  chosen 
such  that  microplasticity  in  cooling  (and  subsequent  unrecovered  strain  in  heating)  took  place 
when  the  stress  was  as  small  as  10  MPa.  The  value  of  the  microplasticity  module  hm  was 
varied.  Fig.7  presents  the  calculated  by  this  model  (for  hm=2  MPa)  diagram  of  strain 
evolution  at  thermal  cycles.  By  adjusting  of  module  hm  one  can  obtain  different  rates  of 
thermocyclic  strain  production  (fig.8).  Simulations  have  also  shown  a  decrease  of  Ms 
temperature. 


Fig.  7.  Strain  dependence  on 
temperature  at  thermal  cycles  in  the 
range  220  -  620  K  under  a  constant 
stress  of  50  MPa  (calculation). 


Fig.8.  Maximum  strain 
accumulated  by  a  model  sample  in  N 
thermal  cycles  under  stress  50  MPa  for 
different  values  of  module  hm :  2  (1),  5 
(2),  7  (3),  65  (4)  MPa. 
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Appearance  of  microplastic  strain  gives  rise  to  some  internal  stress  which  is  obviously 
the  cause  of  the  experimentally  observed  [21]  partial  strain  recovery  at  subsequent  thermal 
cycles  in  an  unloaded  specimen.  Some  confirmation  to  this  statement  is  supplied  by  the  result 
of  modelling  shown  in  fig.9.  Being  very  intensive  for  some  first  cycles,  the  recovery  in  a 
cycle  then  decreases.  In  previous  works  this  effect  was  not  possible  to  describe  by  means  of 
the  structure-analytical  model. 


Fig.9.  Strain  dependence  on 
temperature  at  thermal  cycles  without 
stress  after  preliminary  treatment  shown 
in  fig.  7. 


In  conclusion  we  may  state  that  a  satisfactory  simulation  of  the  cyclic  mechanical 
behavior  of  alloys  with  fccohcp  transformations  is  possible  only  with  an  account  of  the 
accommodation  of  martensite  and  microplastic  deformation. 

Modelling  of  active  plastic  deformation  in  Titanium-nickel  shape  memory  alloy  [4*] 

In  many  applications  there  may  occur  loading  which  causes  an  active  plastic 
deformation.  Besides,  plastic  deformation  influences  the  quantitative  characteristics  and 
kinetics  of  martensitic  non-elasticity  giving  way  to  the  formation  of  such  specific  effects  as 
the  two-way  shape  memory  and  training  of  the  material  at  cyclic  temperature  variations 
[16,21,22].  Preliminary  thermomechanical  treatment  allows  to  change  purposefully  the 
functional-mechanical  properties  of  SMA.  That  is  why  big  attention  is  paid  to  the 
development  of  models  allowing  to  calculate  plastic  deformation  and  predict  its  influence. 
Several  works  are  devoted  to  the  calculation  of  the  dislocation  plasticity  arising  at  the  growth 
of  martensite  crystals  inside  the  austenitic  matrix.  The  principles  of  the  calculation  of  the 
plastic  deformation  due  to  dislocation  slip  have  been  formulated  in  the  book  [9].  In  the 
present  work  it  is  developed  a  model  of  the  plastic  deformation  of  a  polycrystalline  TiNi 
specimen  using  the  general  ideas  of  this  book. 

The  mathematical  object  is  an  austenitic  polycrystal  consisting  of  grains  numbered  1, 
2,...,  co,...,  Q,  each  grain  having  several  slip  systems.  Any  slip  plane  may  be  referred  to  one 
of  1,  2,...,  m,...,  M  types  and  1,  2,...,  k,...,  Km  crystallographically  equivalent  planes  belong 
to  each  type.  These  planes  may  be  obtained  from  one  of  them  by  rotations  belonging  to  the 
point  symmetry  group  of  cubic  lattice.  & 

Athermal  plastic  deformation  ef ' “  in  a  grain  is  the  sum  of  the  deformations  on  each  of  the 
slip  planes  in  this  grain: 


.  M  Km 
m- 1 k~ 1 


(1) 


where  saJm’k)  is  the  referred  to  the  laboratory  basis  athermal  deformation  produced  by  slip 
on  the  £-th  plane  of  the  m-th  type. 

Let  cSr  be  the  effective  stress  applied  to  a  grain.  Then  the  stress  acting  on  the  plane  (m,  k) 
has  shear  components  which  may  be  written  as: 


t0M)  _  Am,k)  J  (m,k)  gr 
T31  ~Ap3  Aq\  Gpq-. 


P( m’k )  _  j(m,k)  j(m,k)  gr 


32 


1 

xp3  Aql 


pq-. 


(2) 


where  Apq  is  the  rotation  matrix  transferring  the  crystallographic  basis  of  the  grain  a  into 
the  crystallographic  basis  of  the  (m,  k)  plane.  We  assume  that  dislocation  slip  starts  when  the 
intensity  of  the  tangential  stress  on  the  slip  plane 


t?" = H^)2) 


^ s(m,k ) 


reaches  its  critical  value  T* (m,k) 

Due  to  the  tangential  stress 

<'*’  =  +«„8l„)  +  T<f»(8l28„3  +  8„8„2) 

a  shear  strain  is  produces,  having  the  components 


(3) 

(4) 

(5) 


where  t<m’k) 
(m,  k). 


na  ( m,k )  _  ( m,k )  _  1  -r(m,k)  X  32 

Psi  2  p  Tx  ’  Ps2  ~  2 Fp  T  ’  (6) 


a  ( m,k)na  ( m,k ) 
pq  P  pq 


is  the  intensity  of  the  shear  strain  rate  tensor  on  the  plane 


We  assume  that  the  flow  stress  Ts{m'k)  is  the  sum  of  an  equilibrium  value  x‘(m)e?,  which  is 
the  same  for  all  planes  of  the  given  m-th  type,  and  additions  responsible  for  the  deformation 
xs{m'k'>  def  and  latent  hardening: 


T*  (m,k)  _  Ti  ( m,k )  eq  +  ^  j  (m,k)  def  +  q  ^  ^  j  (m,*)  def  ^ 

m,k 

where  Ct  is  an  empirical  constant  and  the  last  term  is  due  to  the  latent  hardening  connected 
with  the  possibility  of  simultaneous  plastic  deformation  on  several  slip  planes.  This  formula 
does  not  take  into  account  the  strain  rate  dependence  of  the  flow  stress. 

The  component  xs(m’k)eq  we  shall  write  as  follows: 


^s(m,k)  eq 


T("0 

lF 


k  (m)(T,t^k))T, 


(8) 


where  K(m) (T ,Ip  A))  is  an  empirical  function  specifying  the  temperature  dependence  and 

may  be  a  constant  or  a  distributed  quantity  with  some  density  V|/(m)  (x(m) ) . 

To  calculate  the  evolution  of  the  deformation  hardening  component  of  the  flow  stress 

'Zs(m,k>  def  we  assume  that  the  hardening  modulus  h(m)  depends  linearly  on  the  amount  of 

shear  on  the  plane  so  that  its  rate  Ts  ^m,k^  def  is  equal  to  h(m)t^m,k) .  Besides  hardening 

there  may  take  place  some  recovery  processes  which  lead  to  the  decrease  of  the  flow  stress 
with  the  rate  which  has  an  Arrenius-type  dependency  on  temperature  and  that  of  power-type 
on  the  flow  stress  itself.  Thus 


where  q(T)  =  r(m)ex p 


is{m’k) def  =  -  q(T)( xs(m’k)  de/)m's ,  (9) 

r  Um) 

k  TJ’  >  "V  are  empirical  constants,  k%  is  the  Boltsmann 


constant. 

From  (4)  we  obtain  that  in  the  plastic  state  it  holds: 


_  Ts  (m)  eq  _|_  (m,k)  def  +  q  (m’k)  def  ^ 

m,k 


Taking  time  derivative  of  (10),  using  (8)  and  (9)  and  neglecting  with  the  latent  hardening 
addition  to  the  flow  stress  we  obtain  that  in  the  plastic  flow  state 


•/  ,,  8TS(-m)eq  . 


ST 

or  after  some  transformations 


(ID 


p(  =  rt("*) -K(’\T,t^l>)f  +  q(T  f-’ .  (|2j 

For  P(T , f <"'*))  =P(T),  K<”’(7\f  <"'*>)  =  k<”>  =  const  equation (12) becomes: 


where 


P(T)r^k)  +h{m)t{”  ,k)  =  rp(m,k)  __  (m,k) 

("*>*)  _  K(m)T -  q(T){rs ^ m,k ^ de^ 


(13) 


(13a) 


Assuming  moreover  in  (13)  P(T)  —  0,  and  using  Heaviside’s  functions  to  formulate  the 
flow  condition  (4)  and  the  second  flow  condition  of  the  growth  of  the  difference 

7^>(m,k)  _s(m,k)  .  ,  ^ 

t  “To  ,  we  arrive  at  the  system  of  equations 


p(*M0  = 


j 5  (m>k)  def  =  _  q^T^(xs  ^ def  y 


(14) 


which  solution  allows  to  find  the  rate  of  athermal  plastic  deformation  produced  by  the  action 
of  a  given  slip  system.  The  deformation  of  the  whole  grain  we  get  by  taking  the  sum  (1),  and 
the  macroscopic  plastic  deformation  —  by  neutralization  on  all  the  grains  constituting  the 
polycrystal: 


1  “ 

f  =  -YR 


P  gr  a  n  -1 

©CCD  > 


(15) 


where  Ra  is  the  rotation  matrix  transforming  the  laboratory  basis  into  the  crystallographic  one 
of  the  grain  go. 


Results  of  modeling 


The  approach  described  above  has  been  applied  to  model  plastic  deformation  of 
austenitic  titanium-nickel  alloy.  According  to  the  investigations  of  plasticity  properties  of 
titanium  nickelide  monocrystals  [23]  dislocation  slip  occurres  on  the  crystallographic  planes 
of  two  types  namely  {110}  and  {100}.  To  clear  out  the  contribution  of  the  number  of  active 
slip  systems  there  have  been  carried  out  calculations  of  stress-strain  diagrams  of  a  polycrystal 
for  the  cases  when  only  one  or  both  of  these  system  were  active.  Yield  stress  recovery 
processes  have  been  neglected  (  q(T)  =  0  ),  yield  limits  and  hardening  moduli  have  been 
taken  equal  for  both  slip  systems.  Typical  result  is  presented  in  fig.  10. 


a,  MPa 


Fig.10.  Calculated  stress  -  strain  diagrams  of 
a  model  polycrystal,  consisting  of  20  grains 
when  slip  occurred  on  { 100},  { 1 10}  and  both 
of  these  planes. 


Variation  of  the  hardening  modulus  and  yield  stress  allowed  to  approximate  different 
experimental  curves.  It  may  be  concluded  that  polycrystal  plastic  properties  do  not  depend 
strongly  on  the  number  and  orientation  of  cristallographic  slip  planes.  The  situation  is 


different  for  monocrystals.  Experimental  data  presented  in  the  work  [23]  shows  that  both 
yield  limits  and  hardening  moduli  strongly  differ  and  are  much  bigger  for  a  slip  on  {110} 
planes.  Rater  good  agreement  with  the  experiment  the  present  model  shows  whit  the 

following  values  of  parameters:  for  planes  { 100}  X5 eq  =  70  MPa  ,  h(m)  =  800  MPa  and  for 

planes  {110}  320  MPa  ,  h(m)  -  3000  MPa  (fig.  11).  At  the  same  time  there  is  no 

agreement  when  the  tensile  axis  has  an  orientation  <01 1>  probably  because  of  the  anisotropy 
of  the  flow  stress  for  different  directions  in  the  slip  plane  which  was  not  accounted  in  this 
work. 


cr.MPa 


Fig.  11  .Experimental  (1,2)  [23]  and  calculated 
O’, 2’)  tensile  diagrams  of  titanium  nickelide 
monocrystals  with  the  tensile  axis  oriented 
along  <100>  (curves  1,1’)  and  <110> 
(curves  2,2’) 


Plastic  deformation  when  produced  in  austenitic  state  of  an  SMA  such  as  TiNi  can  be  the 
cause  of  the  reversible  shape  change  when  a  specimen  is  thermocycled  across  the  interval  of 
the  phase  transformations.  This  phenomenon  may  be  explained  by  the  action  of  special 
intergranular  stresses  which  appear  in  a  polycrystal  because  of  the  plastic  deformation 
anisotropy.  An  estimation  calculation  of  the  mean  values  of  these  stresses  added  to  the 
approach  described  above  shows  a  good  qualitative  agreement  with  the  experiment  (fig.  12.). 
These  results  are  soon  to  be  published. 
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Fig.  12  Strain  variation  at  thermocycling  of  a  TiNi  specimen  after  an  active  plastic 
deformation  at  600K: 


a  -  calculation;  b  -  fragment  of  fig.  12a  belonging  under  the  transformation  temperature 
interval  compared  with  experimental  data. 


These  results  show  that  the  structure-analytical  model  of  strength  outlined  by  V.A.Likhachev 
[9]  is  a  good  tool  to  describe  the  plastic  behaviour  of  specific  crystalline  materials  provided 
the  mechanism  of  plasticity  (slip  planes,  critical  stresses,  hardening  moduli)  are  known. 
Particularly,  the  approach  allows  to  model  the  deformation  of  shape  memory  alloys  in 
different  regimes  of  loading  from  the  unique  positions. 

Experimental  study  of  the  influence  of  plastic  strain  on  the  shape  memory 

effect  of  NiTi  [2*] 

The  material  utilized  in  this  study  was  near  equiatomic  NiTi  wires  in  both  the  as-received 
condition  and  after  annealing  at  850°C  for  30  min.  A  differential  scanning  calorimeter  (DSC) 
was  used  to  determine  the  transformation  temperatures  and  latent  heat  for  the  initial  condition 
of  the  specimen.  The  mechanical  loading  portion  of  the  test  consisted  of  an  isothermal  tensile 
loading,  at  a  temperature  above  austenite  finish,  that  imparts  a  measurable  plastic  deformation 
into  the  specimen,  as  seen  in  the  stress  vs.  strain  curve  of  Fig.  13.  After  loading,  the  specimen 
was  heated  in  a  stress  free  condition  to  separate  the  plastic  strain  from  any  stress-induced 
martensite  that  may  have  remained  from  the  mechanical  loading.  This  is  followed  by  a 
thermal  cycle,  under  a  constant  applied  stress,  resulting  in  a  strain  vs.  temperature 
relationship,  Fig.  14,  from  which  the  phase  transformation  characteristics  of  start  and  finish 
transformation  temperatures  and  transformation  strain  can  be  determined.  The  sequence  of  a 
mechanical  loading  followed  by  thermal  loading  was  then  repeated  to  observe  a  change  in  the 
phase  transformation  characteristics  due  to  the  increased  level  of  plastic  strain.  This  entire 
process  was  then  repeated  with  the  isothermal  mechanical  loading  at  a  temperature  below 
martensite  finish.  This  was  performed  to  distinguish  between  the  effects  of  plasticity 
developed  during  phase  transformation  and  plasticity  developed  during  martensitic 
detwinning. 
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Fig.  13  Stress  vs.  Strain  for  as-received  SMA  wire  at  T>  Af 


Temperature,  C 

Fig.  14  Strain  vs.  Temperature  at  25  MPa  after  9.0  %  total  plastic  strain 


Prev.  pi.  def.,  % 

Fig.  15  Dependence  of  phase  transformation  temperatures 
due  to  the  previous  plastic  deformation. 


Fig.  16  The  transformation  temperature  start  and  finish 
differential  due  to  the  previous  plastic  deformation. 


The  following  results  represent  a  complete  series  of  ten  mechanical  and  subsequent 
thermal  loadings  at  a  temperature  above  austenite  finish,  Af.  Fig.  15  demonstrates  influence  of 
plastic  deformation  on  the  characteristic  temperatures  of  the  phase  transformation,  measured 
from  tests  similar  to  Fig.  14  for  an  as-received  SMA  wire.  This  data  demonstrates  that  the 
plastic  strain  significantly  effects  the  martensitic  transformation  temperatures  and  has  a  lesser 
impact  on  the  austenitic  transformation  temperatures.  Fig.  16  quantifies  the  impact  of  the 
plastic  strain  on  the  transformation  temperatures  showing  that  differential  between  the 
martensitic  start  and  finish  temperatures  grow  at  a  rate  more  than  double  that  of  the  austenitic 
start  and  finish  temperatures. 
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Modelling  of  martensite  accomodation  effect  on  mechanical  behavior  of  shape  memory  alloys. 
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An  approach  has  been  presented  to  account  for  micro  plastic  deformation  and  stress  produced  by  accommodation  of 
martensite.  This  has  made  possible  to  describe  such  phenomena  as  incomplete  recovery  of  strain,  strain  accumulation  at 
thermocycling  and  repeated  two-way  shape  memory  effect.  Results  of  modelling  are  in  good  qualitative  agreement  with 
experimental  data. 


Introduction 

At  present  the  structure-analytical  theory  of  strength  (Likhachev  and  Malinin,  1993)  is  successfully  used  to  describe  the 
mechanical  behavior  of  alloys  possessing  the  martensitic  channel  of  non-elasticity.  This  theory  postulates  the  deformation 
law  on  the  microlevel,  and  the  macro  strain  is  calculated  by  neutralization  of  all  micro  strains  related  to  the  volumes  with 
different  crystallographic  orientations  and  parameters  of  the  transformation.  Calculations  have  shown  that  this  approach 
well  describes  the  transformation  plasticity,  shape  memory,  pseudoelasticity  and  some  other  properties.  However  there  has 
not  been  achieved  a  satisfactory  description  of  shape  memory  alloys  (SMA)  behavior  at  cyclic  temperature  variations  (the 
decrease  of  strain  underrecovery  with  the  number  of  cycles)  and  other  properties,  due  to  micro  plastic  strains  produced  at 
the  growth  of  martensitic  crystals,  among  them  the  repeated  shape  memory  effect  (Likhachev  et  al.,  1987).  The  present 
work  is  an  attempt  to  introduce  into  account  this  type  of  micro  plasticity. 


Model 

The  averaging  procedure  as  in  the  work  of  Erglis  et  al.,  (1995)  has  been  organized  in  a  way  to  reflect  the  grain  structure  of 
the  material.  We  have  assumed  that  in  each  of  1,2, ,.,/2grains  of  the  polycrystalline  object  N  crystallographically 
equivalent  martensitic  variants  can  appear.  If  the  matrix  D  of  the  deformation  due  to  the  transformation  by  one  of  the  vari¬ 
ants  is  known,  the  corresponding  matrices  Dn  for  other  variants  can  be  found  by  rotating  D  by  the  operators  Pn  belonging 
to  the  cube  symmetry  group: 

Dn  =  PnDPn-l,  n  =  N .  (1) 


The  volume  fraction  of  martensite  in  a  grain  is  equal  to 


<D 


gr 


1  N 


(2) 


For  simplicity  OSr,  &rb  sSr,  etc  here  and  further  are  written  instead  of  OSr(co),  On(ox),  &r(co),  etc.  Here  the  quantity 
0/N)6n  is  the  volume  fraction  of  the  n- th  variant  of  martensite  in  the  considered  grain. 

In  each  grain  the  average  strain  sSr  accumulated  in  the  course  of  the  transformation  related  to  the  crystallographic 
basis  is 


<3> 

The  macroscopic  strain  t!*1  of  the  specimen  due  to  the  transformation  we  obtain  by  reducing  matrices  (co)  to  the  labo¬ 
ratory  basis  and  averaging  on  all  the  grains : 


1  n 


•f*;1 


(4) 


where  R(a>)  is  the  rotation  matrix  transforming  the  laboratory  basis  into  the  basis  of  grain  ox  * 

To  formulate  the  law  of  variation  of  On  we  consider  the  thermodynamical  equilibrium  of  the  two-phase  austenite- 
martensite  system.  We  introduce  the  thermodynamic  force  Fyf  which  is  the  difference  of  specific  Gibbs  energies  of 
martensite  and  austenite.  It  can  be  split  into  thermal  Fyfh  and  mechanical  Fnmec^  parts: 


Fnt  =  Fyltem  +  Fnmech  pjem  =  TJ(q/TJ,  Fnmech  =  r-Dn  (5) 

Here  T0  is  the  temperature  of  the  equilibrium  of  the  two  phases  when  there  is  no  stress,  qo  is  the  latent  heat  (qo<0),  r  is 
the  stress  and  ••  denotes  the  contraction  of  tensors.  Introducing  the  effective  temperature  Tn*  =  T+  (TJq^vDn,  we  get 
the  expression  for  Fn{  comprizing  both  thermal  and  mechanical  components: 


(6) 


Fnt={f*-T,y{q'/Ta). 


It  is  generally  assumed  that  a  moving  interface  overcomes  a  resistance,  due  to  the  energy  of  martensite  crystals  nu- 
cleation  and  other  barriers.  We  shall  refer  to  the  force  responsible  for  the  hysteresis  as  to  the  friction  force  Fffr  and  assume 
that  it  has  a  constant  absolute  value  and  opposes  the  movement  of  a  phase  boundary: 

F/=  l/^lsignfzld^  (7) 

(we  do  not  put  minus  sign  in  the  right  part  of  (7)  because  of  the  form  in  which  equilibrium  equation  (9)  will  be  written 
further).  Besides  Fyfr,  a  resistance  force  Fn  due  to  interphase  stresses  arises  when  a  martensite  crystal  is  growing.  We 
postulate  as  the  first  approximation  that  this  force  is  proportional  to  On  : 

Fn  =  GmOn>  (8) 


where  Gm  is  a  constant.  We  note  that,  reflecting  the  kinetics  of  the  transformation,  equation  (8)  applies  both  to  cases  when 
martensite  grows  and  when  it  shrinks.  At  the  thermodynamic  equilibrium  there  is  the  balance  of  forces: 

Fnt  =  Fn  +  F/F  (9) 


When  law  (8)  is  specified  constants  |  F.fr  |  and  Gm  can  be  determined  from  the  kinetics  of  the  martensitic  transformation 
when  no  stress  is  applied  : 

(a)  at  T*  =  Ms  On  =  0  and  from  (6)  —  (9)  we  get  |  Ffr  \  =  (~qo)'(Af  -Ms)/(Af+  Ms); 

(b)  at  T*  =  Mf  6n  =  1  and  from  (8),  (9)  it  follows  that  Gm  =  {-qo)-2{Ms  Ms)  . 

If  we  accept  for  phase  equilibrium  temperature  an  estimation  Tg=  ( Af+  Ms)/2  then  we  may  rewrite  the  last  formulas  as 
\Ffr\  =  (^oHAf-Ms)/(2To) ,  Gm  =(-^0H  Ms  -Mj)/T0 . 

Equation  (3)  one  can  consider  as  a  decomposition  of  the  strain  s  Sr  relative  to  some  generalized  basis  {Dn}  with  com¬ 
ponents  On  and  which  we  therefore  shall  interprete  as  the  measures  of  strain.  As  the  micro  plastic  strain  smP  caused  by  the 


accommodation  of  martensite  is  connected  with  the  same  martensite  crystals  as  the  phase  strain,  we  postulate  that  smP  can 
be  expressed  in  the  form  similar  to  (3). 


n  =  J 


(10) 


where  dyP  is  the  measure  of  the  micro  plastic  strain  produced  at  the  accommodation  of  crystals,  belonging  to  the  n- th  vari¬ 
ant  of  martensite  in  a  given  grain. 

As  the  accommodation  leads  to  a  partial  relaxation  of  interphase  stresses,  we  assume  that  the  resistance  force  Fn  de¬ 
creases  by  the  quantity  Gm6fP .  So  we  must  change  equation  (8)  for  resistance  force  to 


Fn=Gm{0n-6rp)  (11) 

We  note  that  after  (1 1)  has  been  accepted  instead  of  (8)  formulas  (a)  and  (b)  for  |  F  fr  |  and  Gm  are  no  more  valid  as 
the  kinetics  of  the  transformation  is  different  when  martensite  accomodation  takes  place. 

To  obtain  an  evolution  for  Oyf  we  assume  that  its  variation  dOyf  is  proportional  to  that  of  the  quantity  of  martensite 
dOn  and  it  occurs  when  the  force  Fn  exceeds  some  yield  limit  F^\ 

d6yP=  k  d0n  7/(|  Fn |  -/y’)  H(dFns\gnFn),  (12) 

where  H  is  the  function  of  Heaviside,  k  is  the  coefficient  characterizing  hardening.  If  k  =  0  there  is  no  micro  plasticity,  if  k 
-»  1  there  is  ideal  micro  plasticity  without  hardening  .  As  the  resistance  force  Fn  grows  (decreases)  only  at  the  growth 
(decrease)  of  martensite,  which  in  its  turn  can  occur  only  at  the  growth  (decay)  of  the  thermodynamic  force  Fnf ,  we  have 
sign (dFn)  =  sign (dFyf)  and 

H(dF  ns\%nFn)=  H{dFnl  signF„).  (13) 


Substituting  (1 1)  into  (9),  passing  to  differentials  and  using  (12  )  and  (13  )  we  arrive  at: 

d0n  =  dFy}/  {Gm[  1-  kH(  \  Fn  \  -F^)H{dFnt  sign F„)]} 


(14) 


Next  we  must  take  into  account  that  micro  plastic  strain  is  accompanied  by  the  hardening  of  material.  To  do  this,  as  it  is 
done  in  the  classical  theory  of  plasticity,  we  assume  that  the  yield  limit  grows  at  the  same  rate  as  the  resistance  force  which 
increment  is  equal  to  that  of  the  thermodynamic  force. 

dFJ  =  dFj  H( I  Fn I  -  F„y  )H(dFnt  sign Fn).  (15) 

To  get  the  initial  conditions  we  use  the  fact  that  at  temperatures  above  v4y  and  in  stress-free  state  no  micro  deformations 
have  begun  and  for  all  martensite  variants  we  have  Fny=  F0y=  const ,  OyP  =  On=  0 . 

When  microplasticity  is  taken  into  account,  the  strain  in  a  grain  is  calculated  as 

'!'=4-i>.K+0»)’  <i6> 

™  n-I 

and  the  total  macro  strain  is  delivered  by  averaging  (4).  Relations  given  above  can  be  reduced  to  a  system  of  Q*N  ordinary 
differential  equations,  which  have  been  used  to  model  the  behavior  a  specimen  made  of  an  alloy  of  TiNi  type  containing 
120  grains  of  different  pseudo-random  orientations.  The  following  values  of  material  constants  have  been  used:  Mf 
=  280  K,  Ms  =  300  K,  As  =  340  K,  Af=  360  K,  To  =  330  K,  qo  =  -150  MJ/m3.  The  value  F(jy  =  13  MPa  has  been  chosen  in 
such  a  way  that  microplasticity  could  take  place  at  cooling  of  specimen  stressed  by  30  MPa.  Parameter  k  has  been  varied. 
Calculations  have  been  performed  for  cases  of  pure  shear  and  uniaxial  tension.  Heat  expansion  has  not  been  taken  into  ac¬ 
count. 

Modelling  of  cyclic  loading 

Microplastic  deformation  at  accommodation  of  martensite  strongly  influences  mechanical  behavior  of  shape  memory 
alloys  and  especially  at  cyclic  variations  of  temperature  or  stress.  Because  of  this  we  have  chosen  to  model  strain  accumu¬ 
lation  and  recovery  by  a  specimen  loaded  with  a  constant  stress  of  100  MPa  and  subjected  to  temperature  cycles  through 
the  full  interval  of  the  phase  transformation,  Fig.l  presents  the  kinetics  of  this  process  for  k  =  0.4.  Maximum  strain  satu¬ 
rates,  strain  unrecovered  in  a  cycle  decreases  and  the  hysteresis  loop  stabilizes  in  approximately  20  cycles  (Fig.2). 


k=0.9 


Fig.l  Strain  vs.  temperature  dependence  at 
thermocycling  in  the  full  interval  of  the 
martensitic  transformation  under  the  stress 
of  100  MPa.  Numerals  on  the  plot  indicate 
numbers  of  cycles. 


Fig.2.  Calculated  strain  at  a  temperature  below  Mf  (a)  and  the  unrecovery  of  the  strain  (b)  in  an  N-th  cycle. 

Thermal  cycling  was  carried  out  through  the  full  interval  of  the  martensitic  transformation  under  the  stress  of 
100  MPa  symbols  in  fig.2b  correspond  to  the  values  of  k,  given  in  fig.2a. 

An  important  place  in  the  studies  of  the  work  ability  of  SMAs  is  occupied  by  the  following  thermomechanical  cycle:  cool¬ 
ing  of  a  sample  from  austenite  to  martensite  under  some  load  P,  application  of  an  additional  load  Q ,  heating  to  a  tempera¬ 
ture  higher  than  AfmA  removing  the  additional  load  Q  to  have  the  initial  load  P.  Figures  3, 4  show  the  result  of  such  nu¬ 
merical  experiment  when  P  -  100  MPa,  Q  =  200  MPa.  Stabilization  of  the  loop  in  this  case  is  also  achieved  in  about  20 
cycles.  Strain  unrecovered  in  a  cycle  becomes  small  and  almost  constant.  The  maximum  strain  accumulated  up  to  the  N-th 
cycle  slowly  grows,  these  results  are  in  a  good  agreement  with  the  data  of  the  direct  experiment  (Belyaev  et  al.,1987).  An 
application  of  an  additional  load  Q  at  a  low  temperature  causes  no  unleastic  strain  as  the  stress  produced  by  P+Q  is  insuffi¬ 
cient  for  the  reorientation  of  martensite.  As  the  temperature  rises,  the  reorientation  limit  stress  decreases.  At  some  tem¬ 
perature  lower  that  yf/the  reorientation  of  martensite  starts.  Then  a  typical  “hump”  appears  on  the  strain  recovery  curves  in 
fig-3. 


Fig.  3.  Calculated  strain  versus  tempera¬ 
ture  diagram  at  thermocycles  with  cooling 
under  the  stress  100  MPa  and  heating  un¬ 
der  300  MPa  (torsion).  Numerals  indicate 
numbers  of  cycles. 


Fig.  4.  Strain  unrecovery  (a)  and  maximum  strain  (b)  as  function 
of  the  number  of  thermocycles  depicted  in  fig.  3. 

Modelling  of  repeated  two-way  shape  memory. 

One  of  the  properties  of  SMAs  which  until  now  could  not  be  satisfactorily  described  in  the  scheme  of  the  structure- 
analytical  model  is  the  repeated  two-way  shape  memory  (Likhachev  et  al.,  1987)  formed  by  cooling  under  a  stress  or  by 
active  straining  in  martensite.  The  mechanism  of  this  phenomena  may  lie  in  the  appearance  of  micro  stress  fields  as  the 
result  of  the  incompatible  plastic  strain  accompanying  the  growth  of  martensite  crystals.  Figures  5  and  6  present  the  results 
of  modelling  of  this  effect. 


Fig.  5.  Calculated  strain  versus  temperature  diagram  at  Fig.  6.  Strain  versus  temperature  diagram  of  active 
cooling  and  heating  under  100  MPa  (curve  1),  unloading  deformation  of  a  sample  in  martensite  state  and  un- 
and  a  stress-free  thermocycle  (curve  2).  loading  (curve  1),  heaning  (curve  2)  and  subsequent 

stress-free  thermocycle  (curve3). 

The  model  demonstrates  a  reversible  strain  variation  which  corresponds  to  experimental  data  (Likhachev  et  al.,  1987).  Ac¬ 
cording  to  existing  notions  (Kuzmin,  1980;  Funakubo  et  al.,1987)  plastic  deformation  at  martensite  accommodation  plays 
an  important  role  in  training  of  the  material  by  thermocycling  under  a  load,  fatigue  damages,  formation  of  eigenstresses 
responsible  for  adaptation  of  the  material  and  formation  of  the  repeated  two-way  shape  memory.  The  last  thesis  is  con¬ 
firmed  by  the  fact  that  the  repeated  shape  memory  appears  only  after  a  strong  enough  straining  of  the  object  by  a  load  ex¬ 
ceeding  a  critical  level  when  there  is  no  complete  strain  recovery. 

At  the  same  time  all  the  previous  considerations  were  mostly  pure  qualitative  and  did  not  allow  a  passage  to  engineer¬ 
ing  calculations.  The  proposed  model  sets  a  formal  ground  for  such  notions  and  the  agreement  of  calculated  results  with 
experimental  ones  confirms  the  validity  of  this  approach.  As  the  result  it  became  possible  to  model  thermocyclic  training 
and  repeated  two-way  shape  memory.  Calibration  of  the  material  constants  allows  to  obtain  valid  results  for  specific  alloys. 
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ABSTRACT 

An  account  of  crystallographic  features  of  fcc<-»hcp  transformations  as  well  as  of  plastic  accommodation 
of  martensite  has  allowed  to  describe  accumulation  of  strain  under  a  constant  stress  at  cyclic  variations  of  tem¬ 
perature.  The  model  also  gives  the  recovery  of  this  strain  at  subsequent  thermal  cycles  when  the  stress  is  re¬ 
moved. 

Keywords:  FeMn,  martensitic  transformation,  modelling,  plastic  accommodation,  plasticity,  thermocycles, 
shape  memory  effect,  unelastic  strain. 


1.  INTRODUCTION. 

One  of  the  specific  features  of  fccohcp  transformations  is  the  multi-variance  of  both  direct  and  reverse 
transformations  '• 2.  The  result  of  it  is  that  specimens  loaded  by  a  constant  stress  experience  deformation  when 
they  undergo  either  of  the  transformations.  Thus,  considerable  strain  can  be  accumulated  at  thermal  cycles. 
Large  part  of  this  strain  is  reversible  and  can  be  recovered  at  further  cycles  after  the  stress  is  removed  This  fact 
suggests  that  the  strain  is  due  to  transformations  rather  than  to  dislocation  movement.  Strain  accumulation  at 
thermocycles  can  be  described  by  the  structure-analytical  theory  of  strength  4  which  is  distinguished  by  two 
ideas:  (1)  the  existence  of  limits  of  the  size  of  growing  martensite  crystals  and  (2)  the  possibility  that  the  strain 
of  a  martensite  crystal  at  the  reverse  transformation  is  not,  in  general,  opposite  to  the  strain  which  this  crystal  has 
undergone  at  the  direct  transformation.  The  model  based  on  only  these  two  features  predicts  absolutely  symmet¬ 
rical  behavior  of  a  loaded  sample  in  heating  and  in  cooling.  An  attempt  to  account  for  an  asymmetry  has  been 
made  in  5.  However,  approaches  4' 5  did  not  take  into  account  microplastic  strain  arising  at  the  growth  of  marten¬ 
site  crystals  and  leading  to  a  partial  accommodation  of  martensite  and  relaxation  of  interphase  stresses.  Because 
of  this  drawback  the  stabilization  of  the  deformation  of  the  model  material  occurs  in  3  -  5  cycles  rather  than  in 
several  dozens  as  in  an  actual  experiment.  Strain  recovery  at  thermal  cycles  without  stress  can  not  be  described 
either.  The  way  to  account  for  microplastic  deformation  in  alloys  with  transformations  of  TiNi  type  has  been 
suggested  in6.  In  the  present  work  an  attempt  has  been  made  to  apply  this  methodology  to  alloys  with  fee  o  hep 
transformations.  Results  of  simulation  of  their  mechanical  behavior  at  thermal  cycles  obtained  by  the  mentioned 
models  are  shown. 


2.  MODEL  WITH  MARTENSITE  CRYSTAL  SIZE  LIMITATIONS 

In  consistence  with  the  ideas  (1)  and  (2)  the  equation  determining  the  evolution  of  the  amount  of  martensite 
belonging  to  any  of  the  orientation  variants  of  lattice  deformation  D„  should  be  written  in  the  form 

H«>H0  -0M00- 

s  f  ■* 

s  f  j 

the  first  term  answers  for  the  direct  and  the  second  for  the  reverse  transformation.  Tn*  =  T+  (To  /  qo )  (x:Dn)  is 
the  effective  temperature  for  the  «-th  variant;  Ms,  Mf,  Af  As  are  the  characteristic  temperatures  of  martensitic 
reactions;  H  is  the  Heaviside’s  function  (//(0)=1);  ®M  =(MN)  Z<D„  is  the  volume  fraction  of  martensite  obtained 
by  averaging  on  all  variants;  dot  denotes  time  derivative;  t  is  the  stress  tensor  and  qo  (normally  negative)  is  the 
latent  heat  of  the  transformation.  Multiplicands  //(<X>°-On)  and  tf(<J>n-<DH)  show  that  direct  and  reverse  transfor¬ 
mations  respectively  are  forbidden  when  the  amount  of  martensite  of  the  «-th  variant  is  greater  than  <t°  or 
smaller  than  <I>H  .  To  reflect  the  second  idea  we  put  OH<0  allowing  negative  values  for  <J>n  meaning  that  amount 
of  austenite  recovered  from  martensite  by  deformation  —  D„  exceeds  the  amount  of  martensite  previously  ob- 


tained  by  deformation.  In  this  model  N  is  the  number  of  all  variants  in  all  the  grains.  Macro  strain  is  obtained  by 
averaging: 

e-(l/A0SOnZ)n.  (2) 

Here  integral  should  be  used  for  continuous  distribution  of  variants  orientations.  According  to  6  satisfactory  re¬ 
sults  are  obtained  with  O0=5,  <X>H=  -4.  Besides  in  4  it  has  been  shown  that  a  better  agreement  with  experiments 
can  be  obtained  when  the  statistic  scatter  of  such  parameters  as  the  strain  temperature  hysteresis  width  and  its 
center  are  taken  into  account.  For  modelling  we  took  the  following  values  of  material  parame¬ 
ters:  ®°“5,  Oh=-4,  q0=-6 3  MJ/m\  D3I  =  DI3  =  0,075  (all  the  rest  £>/*=  0),  Mf=  320  K, 
Ms  ~  370  K,  As  =  470  K,  Af=  520  K.  Maximum  deviations  of  the  hysteresis  width  and  center  from  their 
mean  values  were  taken  75  K  and  50  K  respectively.  All  calculations  were  made  for  uniaxial  tension.  Elastic  and 
heat  expansion  strains  were  not  added.  Typical  diagrams  of  strain  accumulation  in  a  specimen  at  thermocycles 
across  the  interval  200  K  -  620  K  are  shown  in  fig,  1 .  Strain  grows  both  at  the  direct  and  reverse  transformations, 
its  amount  accumulated  in  one  cycle  decreases.  Number  of  cycles  to  stabilize  the  strain  was  the  same  for  simula¬ 
tions  shown  in  figs,  la  and  lb,  more  cycles  needed  for  bigger  values  of  stress  (fig.2). 
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Fig.l.  Dependencies  of  strain  in  the  model  object  under  a  dead  load  150  MPa  on  temperature  at 
thermal  cycles  between  200  K  and  600  K  with  initial  state  in  austenite  (a)  or  martensite  (b). 

These  results  qualitatively  agree  with  direct  observations  of  iron-manganese  alloys 2.  In  an  actual  ex¬ 
periments  more  cycles  are  needed  to  reach  stabilization.  Besides,  calculations  have  shown  that  tempera¬ 
ture  Ms  measured  at  the  start  of  straining  decreases  (fig.3),  this  fact  often  explained  entirely  by  the  growth 
of  the  dislocation  density 7. 

Thus  this  approach  has  demonstrated  the  possibility  of  a  satisfactory  simulation  of  deformation  at 
temperature  variations.  We  note  that  only  strain  due  to  transformations  was  taken  into  account. 


Fig.2.  Dependence  of  the  number  of  cycles  Fig.3.  Temperature  of  the  start  of  strain- 

to  strain  stabilization  on  stress  at  thermocycling  ing  versus  number  of  thermal  cycles  under  a 

across  the  interval  of  the  transformation.  constant  stress  150  MPa. 

3.  MODEL  INCLUDING  MICROPLASTIC  DEFORMATION 

At  fcc-»hcp  transformation  one  of  {111  }fcc  planes  becomes  (000 1  )hcp  by  one  of  three  shears  1/6  <  1 12  >  fcc  on 
each  second  plane  making  up  N=  12  variants  of  transformation  strain  tensor  Dn  and  martensite  quantities  . 

At  the  reverse  transformation  each  of  three  shears  1/3  <  1 120  >hcp  restores  the  initial  orientation  of  austenite. 
Thus  all  variants  can  be  divided  into  four  triplets  (zones)  with  parameters 


(3) 


j  3z 

=~  z=1»2,3,4, 

n=3z-2 

characterizing  the  amount  of  martensite  belonging  to  a  zone,  (1/4)  being  the  volume  fraction  of  martensite  of 
this  zone.  Parameters  <&n  themselves  do  not  have  physical  meaning  of  martensite  amounts  as  deformations  Dn 

belonging  to  one  zone  produce  martensite  of  the  same  orientation.  Still  play  the  role  of  measures  of  the 
phase  deformation. 

As  it  was  mentioned  above  an  incorporation  into  the  model  of  microplastic  strain  caused  by  accommodation 
of  martensite  is  necessary  to  reach  a  quantitatively  correct  modelling  of  cyclic  dependencies.  We  postulate  that 
microplastic  strain  can  be  presented  in  a  form  similar  to  the  lattice  transformation  strain  and  introduce  measures 
of  this  strain  so  that  the  total  micro  strain  p  in  a  grain  is  given  by  the  formula 

^  =  t ri>.  +  <W  (4) 

**  n= 1 

To  formulate  the  kinetics  law  for  Q>n  we  introduce  the  thermodynamic  driving  force 

n=(T'-TayqjT,).  (5) 

which  is  equal  to  the  Gibbs  potential  change  at  the  transformation  of  a  unit  volume  of  austenite  into  martensite 
of  the  n-th  variant.  When  transformation  is  in  process  a  moving  interface  experiences  a  resistance  because  of 
energy  barriers  of  martensite  crystals  nucleation  and  other  obstacles.  The  corresponding  force  responsible  for  the 
existence  of  the  hysteresis  we  settle  to  call  the  friction  force  Ffr  and  assume  that  it  has  a  constant  absolute  value 
and  hinders  the  transformation: 


F"  =  -^"IsignCA®.)  (6) 

Another  opposing  force  Fn  is  that  due  to  the  appearance  of  interphase  stresses.  It  grows  with  and 
is  relaxed  by  microplastic  deformation  .  Thus  we  assume: 

Fn  =Gm(®n  -$«)  (7) 

Gm  is  a  material  parameter.  In  terms  of  these  forces  thermodynamic  equilibrium  equation  is 

(8) 

which  we  used  to  determine  the  conditions  allowing  the  transformation. 

Evolution  equation  for  O^P  we  take  in  the  form  analogous  to  that  in  the  classical  theory  of  plasticity  where 
driving  thermodynamic  force  stands  for  the  stress.  Microplastic  flow  occurs  when  force  Fn  excesses  a  yield 
limit  Fny  and  \f„\  continues  to  grow.  Thus,  when  H(\Fn\  -F„y)  H(dFnsignFn)>0  we  postulate 

dOn=(l//tJdF„  (9) 

with  some  module  of  microplasticity.  Combining  (7)  and  (9)  we  obtain  an  evolution  equation  for  when  mi- 
croplastic  flow  is  on: 


G  *f  fa 

m  rtm 

One  must  keep  in  mind  that  microplastic  deformation  may  be  accompanied  by  hardening.  To  account  for  this  we 
assumed  that  the  yield  limit  increment  d F„y  is  proportional  to  d :  d Fny  =  h'md<S>„P  (h'm  is  some  extra 
material  parameter).  Prior  to  any  plastic  deformations  Fny  =  F0y  =  const  for  all  variants  of  martensite. 

In  iron  manganese  alloys  we  have  two  extra  conditions  of  transformation  finish:  direct  transformation  can 
not  occur  to  make  the  volume  fraction  of  martensite  in  a  grain  <&Sr>  1 ;  martensite  fraction  in  a  zone  (1/4)  ®z  can 
not  become  less  than  zero.  Summarizing  all  this  we  arrive  at  a  final  form  of  the  evolution  equation: 


d F  = 
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Gm  +8», 

„  GmSm 


■  » (K I  -  (-dF,'signF„ )  +  H(F>  -  \ Fn  |)]df. 


+  (^T!LH¥-\ W  -  F^F,) ,  (ii) 

Hi  rim  i 

For  numerical  simulations  the  initial  value  of  the  yield  limit  (^=10  MPa)  was  chosen  such  that  micro¬ 
plasticity  in  cooling  (and  subsequent  unrecovered  strain  in  heating)  took  place  when  the  stress  was  as  small  as 
10  MPa.  The  value  of  the  microplasticity  module  hm  was  varied;  h'm  was  taken  equal  to  hm  .  Fig.4  presents  the 
calculated  by  this  model  (for  hm= 2  MPa)  diagram  of  strain  evolution  at  thermal  cycles  as  in  fig. la,  stabilization 


Fig.  4.  Strain  dependence  on  tempera¬ 
ture  at  thermal  cycles  in  the  range  220  - 
620  K  under  a  constant  stress  of  50  MPa 
(calculation). 


Fig.  5.  Maximum  strain  accumulated  by 
a  model  sample  in  N  thermal  cycles  under 
stress  50  MPa  for  different  values  of  module 
hm:  2  (1),  5  (2),  7  (3),  65  (4)  MPa. 


of  the  strain  being  more  gradual.  By  adjusting  of  module  hm  one  can  obtain  different  rates  of  thermocyclic 
strain  production  (fig.5).  The  decrease  of  Ms  observed  at  simulations  is  similar  to  that  in  fig.3 

Appearance  of  microplastic  strain  gives  rise  to  some  of  internal  stress  which  is  obviously  the  cause  of  the 
experimentally  observed  8  partial  strain  recovery  at  subsequent  thermal  cycles  in  an  unloaded  specimen.  Some 
confirmation  to  this  statement  is  supplied  by  the  result  of  modelling  shown  in  fig.6.  Very  intensive  in  some  first 
cycles  the  recovery  in  a  cycle  then  decreases.  In  previous  works  this  effect  was  not  possible  to  describe  by 
means  of  the  structure-analytical  model. 


Fig.  6.  Strain  dependence  on  temperature  at 
thermal  cycles  without  stress  after  preliminary 
treatment  shown  in  fig.  4. 


In  conclusion  we  may  state  that  a  satisfactory  simulation  of  the  cyclic  mechanical  behavior  of  alloys  with 

fccohcp  transformations  is  possible  only  with  an  account  of  the  accommodation  of  martensite  and  microplastic 
deformation. 
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ABSTRACT 


Micromechanical  modeling  of  the  plastic  deformation  of  austenitic  titanium-nickel  shape 
memory  alloy  poly-  and  monocrystals  has  been  performed.  The  model  assumes  that  the  material 
consists  of  grains  in  which  specific  to  this  alloy  slip  systems  become  active  when  the  shear  stress 
reaches  a  critical  value.  The  plastic  flow  goes  on  at  stress  the  increase  of  which  is  determined  by 
a  hardening  modulus.  To  find  the  parameters  of  the  model  existing  experimental  data  on  the 
plastic  properties  of  monocrystals  have  been  used.  Calculated  deformation  diagrams  of 
monocrystals  with  different  orientations  of  the  tensile  axis  show  a  good  agreement  with 
experiment. 


INTRODUCTION 


The  behavior  of  shape  memory  alloys  (SMAs)  depends  significantly  on  the  internal 
structure  of  the  material.  In  particular,  plastic  deformation  shows  influence  on  the  quantitative 
characteristics  and  kinetics  of  martensitic  non-elasticity,  gives  way  to  the  formation  of  such 
specific  effects  as  the  two-way  shape  memory  and  training  of  the  material  at  cyclic  temperature 
variations  (Likhachev,  Kuzmin  and  Kamentzeva,  1987;  Belyaev,  Kuzmin,  Likhachev  and 
Kovalyev,1987;  Airoldi,  Ranucci,  Riva  and  Sciacca,1995).  Preliminary  thermomechanical 
treatment  allows  to  change  purposefully  the  functional-mechanical  properties  of  SMAs.  That  is 
why  big  attention  is  paid  to  the  development  of  models  allowing  to  calculate  plastic  deformation 
and  predict  its  influence.  Several  works  are  devoted  to  the  calculation  of  the  dislocation  plasticity 
arising  at  the  growth  of  martensite  crystals  inside  the  austenitic  matrix.  It  has  been  shown  that 
accomodation  processes  are  connected  with  such  phenomena  as  the  decrease  of  the  phase  yield 
limit  at  cyclic  loading  (Patoor,  Siredey,  Eberhardt  and  Berveillerl995),  accumulation  of 
irreversible  strain  at  thermal  cycling  under  a  load  or  without  it  (Bo,  Lagoudas,  1998;  Evard, 
Volkov,  1999).  The  main  principles  of  the  calculation  of  the  plastic  deformation  due  to 
dislocation  slip  have  been  formulated  in  the  book  (Likhachev,  Malinin,  1993)..  The  present  work 
is  an  attempt  using  the  general  ideas  of  this  book  to  develop  a  model  of  the  plastic  deformation 
of  polycrystalline  TiNi  specimen. 


THEORETICAL  MODEL 


The  mathematical  object  is  an  austenitic  polycrystal  consisting  of  grains  numbered  1, 
2,...,  to,...,  Q,  each  grain  having  several  slip  systems.  Any  slip  plane  may  be  referred  to  one  of  1, 
2,...,  m,...,  M  types  and  1,  2,...,  £,...,  Km  crystallographically  equivalent  planes  belong  to  each 
type.  These  planes  may  be  obtained  from  one  of  them  by  rotations  belonging  to  the  point 
symmetry  group  of  cubic  lattice. 


Athermal  plastic  deformation  efa  in  a  grain  is  the  sum  of  the  deformations  on  each  of 
the  slip  planes  in  this  grain: 


M  Km 

Br=zz<“’‘) 

m=l  k- 1 


(i) 


where  is  the  referred  to  the  laboratory  basis  athermal  deformation  produced  by  slip  on 

the  £-th  plane  of  the  m-th  type. 


Let  ct Sr  be  the  effective  stress  applied  to  a  grain.  Then  the  stress  acting  on  the  plane 
(m,  k)  has  shear  components  which  may  be  written  as: 


x<7-‘)  =  ,  *$"•*>  =  A(’MA<r-‘>cZ 
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where  Apq  is  the  rotation  matrix  transferring  the  crystallographic  basis  of  the  grain  e>  into  the 
crystallographic  basis  of  the  (m,  k)  plane.  We  assume  that  dislocation  slip  starts  when  the 
intensity  of  the  tangential  stress  on  the  slip  plane 


reaches  its  critical  value  x 


Due  to  the  tangential  stress 
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a  shear  strain  is  produces,  having  the  components 
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where  fp  ’  }  is  the  intensity  of  the  shear  strain  rate  tensor  on  the  plane 


(m,  k). 


We  assume  that  the  flow  stress  xs(m,k)  is  the  sum  of  an  equilibrium  value  x 


s(m)eq 
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which  is  the  same  for  all  planes  of  the  given  m- th  type,  and  additions  responsible  for  the 
deformation  xs{m'k)  def  and  latent  hardening: 

^  j  (m,k)  _^s  ( m,k )  eq  (m,k)  def  +  q  ^  X^  (m,i)  ^  ^ 

m,k 

where  C  T  is  an  empirical  constant  and  the  last  term  is  due  to  the  latent  hardening  connected  with 

the  possibility  of  simultaneous  plastic  deformation  on  several  slip  planes.  This  formula  does  not 
take  into  account  the  strain  rate  dependence  of  the  flow  stress. 

The  component  xs(m'k)eq  we  shall  write  as  follows: 

x' (-•*)*  =  x?) -K(,,(r,fpW))T,  (8) 

where  Kw(r,rpw))  is  an  empirical  function  specifying  the  temperature  dependence  and 

(  m\ 

X/r  may  be  a  constant  or  a  distributed  quantity  with  some  density  vj/m)(x(m)) . 

To  calculate  the  evolution  of  the  deformation  hardening  component  of  the  flow  stress 
xs(m,k)def  we  assume  that  the  hardening  modulus  h(m)  depends  linearly  on  the  amount  of 

shear  on  the  plane  so  that  its  rate  Ts  (m,k)  def  is  equal  to  h(m)t^m'k) .  Besides  hardening  there 


may  take  place  some  recovery  processes  which  lead  to  the  decrease  of  the  flow  stress  with  the 
rate  which  has  an  Arrenius-type  dependency  on  temperature  and  that  of  power-type  on  the  flow 
stress  itself.  Thus 
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where  q(T)  =  r,")  exp 
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;  U m  ,  «t,  are  empirical  constants,  k%  is  the  Boltsmann 


constant. 

From  (4)  we  obtain  that  in  the  plastic  state  it  holds: 

_  ^s(m)  eq  _j_  ^  j  (m,k)  def  ^  ( m,k )  def 

m,k 

Taking  time  derivative  of  (10),  using  (8)  and  (9)  and  neglecting  the  time  derivative  of  the  latent 
hardening  addition  to  the  flow  stress  we  obtain  that  in  the  plastic  flow  state 
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or  after  some  transformations 

P(T,t^k))f^m’k)  +h{m)t^k)  =Ttm’k) -Yi(m\Tj^k))f  +  q{T){xs{m'k)def)m^  .  (12) 
For  P(T,t^m,k))  =  P(T),  K{m)(T>t^m'k))  =  K{m)  =const  equation (12) becomes: 

P(7’)fp(m'*)  +  h(m)t^k) 

where 
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Assuming  moreover  in  (13)  P(T)  =  0,  and  using  Heaviside’s  functions  to  formulate  the  flow 

condition  (4)  and  the  second  flow  condition  of  the  growth  of  the  difference  T^m,k^  —  xs0^m,k\  we 
arrive  to  the  system  of  equations 

=(h{m)y'[t}m’k)  -Ts0(m’k) ] 

(m,k)  def  _  ^(m)p(m,*)  _  q(Y)(xs  («»*)  def 

which  solution  allows  to  find  the  rate  of  athermal  plastic  deformation  produced  by  the  action  of  a 
given  slip  system.  The  deformation  of  the  whole  grain  we  get  by  taking  the  sum  (1),  and  the 


macroscopic  plastic  deformation  —  by  neutralization  on  all  the  grains  constituting  the 
polycrystal: 
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where  Ra  is  the  rotation  matrix  transforming  the  laboratory  basis  into  the  crystallographic  one  of 
the  grain  co. 


RESULTS  OF  MODELING 


The  approach  described  above  has  been  applied  to  model  plastic  deformation  of  austenitic 
titanium-nickel  alloy.  According  to  the  investigations  of  plasticity  properties  of  titanium 
nickelide  monocrystals  (Chumlyakov,  Surikova  and  Korotaev,1996)  dislocation  slip  occurres  on 
the  crystallographic  planes  of  two  types  namely  {110}  and  {100}.  To  clear  out  the  contribution 
of  either  of  these  slip  systems  there  have  been  carried  out  calculations  of  stress-strain 
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Fig.l.  Calculated  stress  -  strain  diagrams  of  a 
model  polycrystal,  consisting  of  20  grains 
when  slip  occurred  on  {100},  {110}  and  both 
of  these  planes. 


diagrams  of  a  polycrystal  for  the  cases  when  only  one  or  both  of  these  system  were  active.  Yield 
stress  recovery  processes  have  been  neglected  (  q(T)  =  0  ),  yield  limits  and  hardening  moduli 
have  been  taken  equal  for  both  slip  systems.  Typical  result  is  presented  in  fig.l.  Variation  of  the 
hardening  modulus  and  yield  stress  allowed  to  approximate  different  experimental  curves.  It  may 
be  concluded  that  polycrystal  plastic  properties  do  not  depend 


a,  MPa 


Fig.2  .Experimental  (1,2)  (Chumlyakov, 
Surikova  and  Korotaev,  1996)  and  calculated 
(l’,2’)  tensile  diagrams  of  titanium  nickelide 
monocrystals  with  the  tensile  axis  oriented 
along  <100>  (curves  1,1’)  and  <110> 
(curves  2,2’) 

strongly  on  the  number  and  orientation  of  cristallographic  slip  planes.  The  situation  is  different 
for  monocrystals.  Experimental  data  presented  in  the  work  (Chumlyakov,  Surikova  and 
Korotaev,  1996)  shows  that  both  yield  limits  and  hardening  moduli  strongly  differ  and  are  much 
bigger  for  a  slip  on  {110}  planes.  Rater  good  agreement  with  the  experiment  the  present  model 

shows  whit  the  following  values  of  parameters:  for  planes  {100}  Tseg  =  70  MPa  , 

h(m)  -  800  MPa  and  for  planes  {110}  Tseg  =  320  MPa  ,  h(m)  =  3000  MPa  (fig.2).  At  the  same 
time  there  is  no  agreement  when  the  tensile  axis  has  an  orientation  <01 1>  probably  because  of 
the  anisotropy  of  the  flow  stress  for  different  directions  in  the  slip  plane  which  was  not 
accounted  in  this  work.  . 

Making  the  resume  of  this  work  we  may  note  that  the  structure-analytical  model  of 
strength  outlined  by  V.A.Likhachev  (Likhachev,  Malinin,  1993)  is  a  good  tool  to  describe  the 
plastic  behavior  of  specific  crystalline  materials  provided  the  mechanism  of  plasticity  (slip 
planes,  critical  stresses,  hardening  moduli)  are  known.  Particularly,  the  approach  allows  to  model 
the  deformation  of  shape  memory  alloys  in  different  regimes  of  loading  from  the  unique 
positions. 
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ABSTRACT:  Fully  annealed  wires  of  NiTi  are  loaded  mechanically  under 
isothermal  conditions  until  plastic  deformation  occurs.  After  unloading,  a 
temperature  cycle  under  a  constant  stress  is  applied  to  determine  a  relationship 
between  the  applied  plastic  strain  and  the  phase  transformation  characteristics. 
Two  sets  of  isothermal  mechanical  loadings  are  performed,  with  the  first  at  a 
temperature  below  the  martensitic  finish  temperature,  and  the  second  at  a 
temperature  above  the  austenitic  finish  temperature. 

INTRODUCTION:  Unlike  common  ductile  metals,  plastic  deformation  in  shape 
memory  alloys  can  be  created  by  phase  transformations  and  can  occur  at  relatively 
low  stress  levels.  Along  with  the  creation  of  the  plastic  deformation,  other 
material  properties  under  go  significant  changes  with  the  accumulation  of  plastic 
strain.  The  austenite/martensite  phase  transformation  characteristics  of  shape 
memory  alloys  have  been  shown  in  previous  studies  (Hebda  and  White  [1995], 
Lim  and  McDowell,  [1994],  Liu  and  McCormick,  [1990],  McCormick  and  Liu’ 
[1994],  Perkins  and  Bobowiec  [1986],  Perkins  and  Sponholz  [1984])  to  be  related 
to  the  loading  history  of  the  alloy.  Both  thermal  cycles,  mechanical  cycles  and 
cold  working  have  shown  an  effect  on  the  phase  transformation  thermomechanical 
properties  such  as  start  and  finish  temperatures,  transformation  strain  and  latent 
heat  of  transformation.  The  growth  of  plastic  deformation  has  been  successfully 
modeled  and  predicted  by  Bo  and  Lagoudas  [1998]  during  thermally  induced 
phase  transformations  under  constant  applied  stress.  To  further  develop  and  verify 
this  model,  experiments  are  conducted  to  study  the  influence  of  mechanically 
induced  plastic  deformation  on  the  characteristics  of  the  austenite/martensite 
phase  transformation. 


PROCEDURES,  RESULTS  AND  DISCUSSION:  The  material  utilized  in  this 
study  was  near  equi-atomic  NiTi  wires  in  both  the  as-received  condition  and  after 
annealing  at  850°C  for  30  min.  A  differential  scanning  calorimeter  (DSC)  was 
used  to  determine  the  transformation  temperatures  and  latent  heat  for  the  initial 
condition  of  the  specimen.  The  mechanical  loading  portion  of  the  test  consisted  of 
an  isothermal  tensile  loading,  at  a  temperature  above  austenite  finish,  that  imparts 
a  measurable  plastic  deformation  into  the  specimen,  as  seen  in  the  stress  vs.  strain 
curve  of  Fig.  1.  After  loading,  the  specimen  was  heated  in  a  stress  free  condition 
to  separate  the  plastic  strain  from  any  stress-induced  martensite  that  may  have 
remained  from  the  mechanical  loading.  This  is  followed  by  a  thermal  cycle,  under 
a  constant  applied  stress,  resulting  in  a  strain  vs.  temperature  relationship,  Fig.  2, 
from  which  the  phase  transformation  characteristics  of  start  and  finish 
transformation  temperatures  and  transformation  strain  can  be  determined.  The 
sequence  of  a  mechanical  loading  followed  by  thermal  loading  was  then  repeated 
to  observe  a  change  in  the  phase  transformation  characteristics  due  to  the 
increased  level  of  plastic  strain.  This  entire  process  was  then  repeated  with  the 
isothermal  mechanical  loading  at  a  temperature  below  martensite  finish.  This  was 
performed  to  distinguish  between  the  effects  of  plasticity  developed  during  phase 
transformation  and  plasticity  developed  during  martensitic  detwinning. 

The  following  results  represent  a  complete  series  of  ten  mechanical  and 
subsequent  thermal  loadings  at  a  temperature  above  austenite  finish,  Af.  Fig.  3 
demonstrates  influence  of  plastic  deformation  on  the  characteristic  temperatures 
of  the  phase  transformation,  measured  from  tests  similar  to  Fig.2  for  an  as- 
received  SMA  wire.  This  data  demonstrates  that  the  plastic  strain  significantly 
effects  the  martensitic  transformation  temperatures  and  has  a  lesser  impact  on  the 
austenitic  transformation  temperatures.  Fig.  4  quantifies  the  impact  of  the  plastic 
strain  on  the  transformation  temperatures  showing  that  differential  between  the 
martensitic  start  and  finish  temperatures  grow  at  a  rate  more  than  double  that  of 
the  austenitic  start  and  finish  temperatures. 


Figure  1  Stress  vs.  Strain  for  as-received  SMA  wire  at  T>  Af 


Temperature,  C  3*  Strain 


vs.  Temperature  at  25  MPa  after  9.0  %  total  plastic  strain 


Figure  3  Dependence  of  phase  transformation  temperatures 
due  to  the  previous  plastic  deformation. 
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Figure  4  The  transformation  temperature  start  and  finish 
differential  due  to  the  previous  plastic  deformation. 

CONCLUSIONS:  In  this  study,  experiments  were  performed  to  identify  the 
effect  of  plastic  deformation  on  the  phase  transformation  properties  of  NiTi.  Test 
results  show  that  the  martensitic  transformation  temperatures  widen  as  the  plastic 
strain  increases  in  the  material.  These  results  are  consistent  to  those  found  for 
plastic  strain  development  under  cyclic  thermal  loading  with  constant  applied 
stress. 
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